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O ' ABSTRACT 

o : 

CN| ■ We briefly review the high energy astro-physical processes that are related to 

$_i ' the production of high energy 7-ray and neutrino signals and are likely to be 

Ohi important for the energy loss of high and ultrahigh energy cosmic rays. We also 

' give examples for neutrino fluxes generated by different astrophysical objects 

00 . and describe the cosmological link provided by cosmogenic neutrinos. 

t-h ; 

1. Introduction 

> ; 

The traditional high energy astrophysics is based purely on electromagnetic pro- 
cesses. Its recent extensions include the models for the production of TeV 7-rays, 
which are tested versus observations and proven to be conceptually correct. Such 

^ 1 7-rays are generated by the inverse Compton effect (ICE), in which accelerated elec- 

trons boost seed photons to the TeV energy range. Both basic types of ICE models: 

^ ■ the synchrotron-self Compton model (SSC) and inverse Compton boosting of external 

photons can provide a good description of the observational data. In the SSC model 
seed photons are generated by synchrotron radiation of the accelerated electrons, 

"£0 \ which then boost them to high energy. 

ICE models describe quite well the double peaked structure of the emission of AGN 
jets. The lower energy peak in the sub-MeV range represents the seed radiation, while 
the TeV peak is boosted by the accelerated electrons. The model parameters are also 
quite consistent with the high variability of the signals observed from active galactic 
nuclei (AGN), where the TeV 7-ray flux is known to double on the timescale of 
minutes. In addition, the SSC model gives a specific relation between the amplitudes 
of the TeV and sub-MeV signals, which is at least qualitatively similar to the observed 
one. 

Neutrino astrophysics, on the other hand, is based on hadronic interactions that 
are not proven to happen in any astrophysical system. Defining the role of hadronic 
processes in the dynamics of powerful astrophysical systems is one of the aims of 
the neutrino astrophysics. There are though no reasons to believe that the observed 
TeV 7-ray signals are necessarily generated in purely electromagnetic processes and 
hadronic interactions have no contributions to them. Following the pioneering work 
in the last couple of decades there are now hadronic models that describe equally 
well the double peaked 7-ray energy spectrum. The lower energy photons are again 



generated by synchrotron radiation of electrons, but the electrons themselves are 
results of hadronic production of mesons and meson decays that feed electromagnetic 
cascades. The fast variability of the sources is more difficult to predict in hadronic 
models. 

We shall discuss both types of processes. Neutrino producing processes include in- 
elastic pp and pj collisions, while the non-neutrino producing processes are the electro- 
magnetic (Bethe-Heitler) pj interactions, 77 interactions, and the inverse Compton 
effect. Any astrophysical model has to include synchrotron radiation, which cannot 
be avoided in the relatively high magnetic fields that are needed to accelerate charged 
particles. Synchrotron radiation is the production process for sub-MeV photons, from 
X-rays down to radio waves. 

We will not discuss particle acceleration. In principle, both protons and electrons 
should be accelerated to the same Lorentz factors in all astrophysical systems. Ac- 
celeration processes do not care about the particle charge. Injection, however, could 
be different. The injection of protons is better understood and described in detailed 
shock acceleration studies. 

The same neutrino producing (and not producing) processes are important for 
the energy loss and propagational effects on ultra high energy cosmic rays (UHECR) 
in the microwave background (MBR) and in other radiation fields. Because the 
energy spectra of seed fields are quite different in different objects we will discuss the 
interaction lengths in the MBR, which provides the only standard photon field energy 
spectrum and allows for a reasonable comparison of different processes. 

2. Neutrino producing processes 

High energy astrophysical neutrinos are generated in the decay chains of mesons, 
such as 7T ± — ► ^ + v^iv^), ^ — > + v e (i> e ). Other mesons are also involved 

to certain degree and all decay chains are very well known. What is not known is 
how the mesons are produced. 

Inelastic pp collisions are no doubt the best studied process in high energy physics. 
The proton interaction length in Hydrogen is 51 g/cm 2 , shorter than the radiation 
length of 61 g/cm 2 . For the average nucleon density in the Galaxy (1 cm~ 3 ) this 
converts to a distance of 3xl0 25 cm, i.e. 10 Mpc, almost three orders of magnitude 
higher than the linear size of our Galaxy. If the Hydrogen target were a dense molec- 
ular cloud with density of 300 cm" 3 one interaction length would coincide with the 
diameter of the Galaxy. These numbers refer, of course, to the proton pathlength, 
and not directly to linear dimensions. Protons may be contained in magnetic fields, 
which will increase their interaction probability. 

Because of the small average density of matter in astrophysics there are only a 
few objects that can present targets for pp interactions. These are: 



• stars 



• accretion discs, that can contain column densities of 50 g/cm 2 close to the 
compact object, and 

• dense molecular clouds, compressed, for example, by expanding supernova rem- 
nants. The density of such clouds can reach 1,000 cm" 3 . 

The disadvantage of the pp collisions as important astrophysical process is in the 
rarity of the needed target material. The main advantage is the very low interaction 
threshold - a Lorentz factor of 1.3 for the production of the A resonance. 
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Figure 1: Energy loss length and protons in the microwave background. 



Proton photoproduction is on the opposite end of interaction properties. Its main 
advantage is the availability of photon targets in all astrophysical systems, and even 
in extragalactic space. Only MBR has a number density in excess of 400 cm -3 which 
fully compensates for the cross section ratio u pl /a pp which is of order 0.01. The 
main disadvantage of the photoproduction is its high energy threshold. The center 
of mass energy E p e(l — cos6) should exceed the sum of the proton and pion masses 
squared. Since the photon background fields are often sub-eV this turns out to be a 
hard requirement. 

Another difference between the two processes is the proton inelasticity coefficient 
K ine i which describes the fractional energy loss of the proton in a collision. In in- 
elastic pp interactions K ine i is about 0.5 and grows logarithmically with energy. In 



photoproduction interactions on the MBR Ki ne i is 0.17 for protons of energy 10 20 eV, 
0.27 at 10 21 eV and only asymptotically approaches 0.5. The important astrophysical 
quantity Xi oss , the energy loss length, equals X/K ine i. At 10 20 eV it is almost a factor 
of six longer than the interaction length. 

Figure ^ shows the energy loss length for protons in the microwave background. 
The dashed line is the interaction length for photoproduction X p h and the thin solid 
lines show the Xi oss for photoproduction and for the Bethe-Heitler pair creation pro- 
cess. The thick solid line is the total Xi oss which also includes the adiabatic energy 
loss length of 4 Gpc for H = 75 km/s/Mpc. 

The dotted line shows the neutron decay length. Only neutrons of energy ex- 
ceeding 4xl0 20 eV undergo photoproduction interactions in the MBR. The neutron 
and proton photoproduction cross sections are almost identical, except at the energy 
threshold of the process. Lower energy neutrons always decay and generate i> e , so one 
has to include neutron decay in the list of neutrino producing processes. 

3. Electromagnetic processes 

3.1. Pair creation by protons 

The line marked BH shows the proton energy loss length in the electromagnetic 
pair creation process p'j — > pe + e~ . This is an interesting process that may have a cru- 
cial importance for understanding the fate of ultrahigh energy cosmic rays (UHECR) 
in the Universe since it creates a feature in their spectrum^. 
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Figure 2: Energy spectra of positrons generated by protons of different energies in the MBR. 
The proton energy loss in this process peaks at about 2xl0 19 eV. The cross section 




of the process monotonically grows with the energy and the shape of the energy loss 
curve is determined by the decreasing proton inelasticity. K ine i falls from 7xl0 -4 at 
10 18 eV to 10~ 5 at 10 21 eV. For this reason the electrons of the created pair have a 
similar energy distributions which is independent from the process cross section. 

Figure El shows the spectra of the positrons generated in the MBR by protons of 
energy 10 18 , 10 19 , and 10 20 eV protons. Electron spectra are the same and the total 
energy loss is twice the integral of the histograms shown in Fig. El 

The hardest electrons are generated by 10 19 eV protons. For higher proton energies 
the electron spectrum is wider, but peaks at almost the same energy as at lower proton 
energy. The growing cross section can be visually detected by the increasing area of 
the histograms with the proton energy. 

3.2. Inverse Compton scattering 

Inverse Compton scattering is described by the same expressions as the Compton 
effect. High electrons interact with the seed photons and boost them to higher energy. 
The relation between the energies of the two particles is 

k = ^{l- (3cos0) 
m e c 2 

At low CM energies (eE e <C (m e c 2 ) 2 ) the process is in the Thomson regime which 
is characterized by the Thomson cross section cr^ = 87rr 2 /3 (665 mb) and a flat 
distribution of the boosted photons. At higher energy there is a transition to the 
Klein-Nishina regime during which the cross section decreases and the boosted pho- 
tons become hard. Although the cross section decreases as s -1 / 2 the interaction length 
may have a slightly different behavior as electrons interact on different parts of the 
seed photon spectrum. Asymptotically the boosted photon energy equals the electron 
energy. 

3.3. Gamma-gamma collisions 

In 7-ray producing environments inverse Compton scattering is always accompa- 
nied by the process 77 — > e + e~ which is the opposite to electron-positron annihi- 
lation. The process cross section peaks at s — 4m 2 . In astrophysical settings the 
resonant behavior is less noticeable and the cross section peak is smoothed by the 
wide energy spectrum of the seed photons. Fig. El shows the interaction lengths for 
inverse Compton scattering and for pair production on the microwave background. 

At low energy one can see the ICE interaction length in the Thomson regime. 
The pair production process is below threshold in this energy range. It reaches its 
minimum interaction length at E 1 = 1.20 xlO 6 GeV. Soon after that both cross sec- 
tions start to decline at the same rate. As far as interactions in extragalactic space 
are concerned, the interaction lengths are modified by interactions on photon fields 
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Figure 3: Interaction lengths for ICE (dots) and pair production by 77 collisions in the MBR. 

different from MBR. At energy lower than 10 15 eV the infrared/optical background 
plays a major role. The observations of TeV 7-rays from distant AGN are often used 
to put lim its on the density of the isotropic IR/O background in the near infrared 
region At energies above 10 19 eV interactions on the radio background are im- 
portant. Since the radio background in the MHz region is not observable, there is a 
large uncertainty in the 7-ray interaction length above 10 19 eV. At still higher energy 
comes in the two pairs production process 77 — > e + e~ e + e~ which limits the 7-ray 
interaction length to 80 Mpc even in the absence of radio background. 

The similar cross sections of these two processes and their small interaction lengths 
support the development of electromagnetic cascades in astrophysical objects and in 
the extragalactic space. Photons interact to produce e + e~ pair. Electrons then scatter 
the seed photons to high energy and the scattered photons generate another pair. 
At high energy, when both ICE and pair production interactions have high energy 
transfer, i.e. the secondary particles have almost the same energy as the interacting 
ones, the cascade development leads only to a slow degradation of the initial 7-ray 
energy. Close to the 77 energy threshold the pair electron distribution becomes more 
uniform and the energy degradation speeds up. 

The major contributor to the acceleration of these electromagnetic cascades is the 
synchrotron radiation that high energy electrons suffer in the presence of magnetic 
fields. 

3.4- Synchrotron radiation 

Synchrotron radiation (magnetic bremsstrahlung) is a very important energy loss 



process for charged particles in the presence of magnetic fields. The energy loss is 
proportional to the electron Lorentz factor squared times the energy density of the 
magnetic field and depends on the pitch angle of the electron with the magnetic field 
lines. For an ensemble of relativistic electrons that are scattered randomly in all 
directions the energy loss averaged over all pitch angles in particle physics units is 

dt V Gauss/ VGeV/ 1 

The characteristic frequency of the radiated photons is their critical frequency 

= 161 x 1013 (d=) 

Expressed as a fraction of the electron energy the critical frequency is proportional to 
the product of the energy and the magnetic field value i> c oc E 2 x B. The number of 
emitted photons peaks at 0.29^ c . The higher the energy, the harder is the spectrum 
of the radiated photons as illustrated in Fig. HJ 
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Figure 4: Synchrotron photon power for 10 14 eV electrons in 1 fiG field (left) and for 10 10 eV 
electrons in 1 nG field. 

The total energy loss of a 10 19 eV electron in 10~ 9 G field is higher than that of a 
10 14 eV electron in fiG field by 10 4 , but the critical frequency of the photons it emits 
is higher by 10 7 . Note that synchrotron radiation brings the energy of the 10 19 eV 
electron directly into the GeV range. 

Since the synchrotron energy loss depends on the square of the particle Lorentz 
factor it is thus inversely proportional (for the same total energy) to the square of 
the particle mass. A proton loses only (m e /m p ) 2 ~ 3 x 10~ 7 times as much energy 




as an electron of the same E tot . The energy loss of muons is 2.37xl0 -5 down from 
the electron one for the same energy. Proton and muon energy losses on synchrotron 
radiation can thus be important only in very strong magnetic fields and respectively 
higher particle energy. 

Figure 03 compares the proton energy loss length to the photon 77 interaction 
length (which can be considered energy loss length in the presence of magnetic field) 
in the microwave background. The radio background^ is included in the calculation 
of the photon interaction length. Because of that the energy loss distances for the 
highest energy photons are still highly uncertain. 
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Figure 5: The energy loss length for protons in MBR is compared to the photon energy loss length. 



Electromagnetic processes lead to faster energy loss for high energy photons and 
electrons than nucleons lose in hadronic interactions. The actual interaction rates, 
however, depend strongly not only on the total energy density of the seed radiation, 
but also on its energy spectrum. Since all particle production processes have energy 
thresholds, and only the seed photons above that threshold count, different environ- 
ments could favor electromagnetic or hadronic processes as demonstrated in® in the 
case of BL Lac objects. 

The immediate conclusion from this brief review is that all physics processes are 
very well known, and the main uncertainties in our estimates come from the lack of 



knowledge on the astrophysical environments in potential neutrino sources. 

4. Resonant neutrino cross sections 

It will be a loss not to mention a couple of resonant neutrino cross sections that 
have never been measured, but have been calculated on a solid physics basis. 




Figure 6: Left hand panel: cross section for the 'Glashow' resonance, shown with one of its extensions. 
Right hand panel: neutrino annihilation, and resonant neutrino-photon interaction. The cross 
section does not depend on the neutrino mass. 

The one shown in the left hand panel of Fig. 03 is the resonant interaction of 
u e with electrons, that was suggested by Glashow^. The idea was developed to 
its current understanding by Berezinsky & Gazizov^. The cross section peaks at 
E De = My V /(2m e ) = 6.4 xlO 6 GeV. The maximum value including all nine W~ 
decay channels is 4.7xl0~ 31 cm 2 . There are numerous extensions the resonant v e e~ 
scattering. We show one of them, E=e~ — > , that enhanced the width of the 
resonance to higher neutrino energy^*. 

The right hand panel shows the neutrino annihilation [yv — > Z ) cross section 
which peaks at energy Z^jra v . It is plotted for m v = 1 eV. This process has been 
used as a basis of the Z-burst scenario for the ultrahigh energy cosmic rays, where 
UHECR are result of Z decay within the GZK sphere of about 50 Mpc The 
cross section for the process wy — > e~W + (right hand panel) does not depend on the 
neutrino mass but it does depend on the lepton mass. In the contemporary universe 
it peaks at neutrino energy about 10 16 GeV and is totally insignificant. It does, 
however, lead to absorption of all high energy neutrinos generated at high redshifts. 



E e = 10 10 GeV 
B = 1 nG 



The influence of this process is felt already at [z > 10) in interactions on the MBR 



5. Examples of predicted neutrino fluxes 

We do not have the ambition to review the models for the production of high 
energy astrophysical neutrinos in this talk. The discussion of the relevant processes 
would be, however, not complete without several examples of predicted fluxes. 

5.1. Source neutrinos 

Figure [7| shows + fluxes predicted for five potential neutrino sources. The 
atmospheric neutrino fluxes within 1° from the source are indicated with a shadowed 
region. The upper edge of that region corresponds to horizontal neutrinos and lower 
one - to vertical neutrinos. 

The curve labeled 1) shows the neutrinos that we expect from the direction of 
the Sun^J. High energy neutrinos are generated by cosmic ray interactions in the 
rarefied solar envelope and some of them propagate through it to us. The energy 
spectrum of these neutrinos is flatter than the one of atmospheric neutrinos because 
mesons decay easier in the tenuous environment. 
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Figure 7: Examples of several different source neutrino fluxes. See text for a brief description of the 
models. 

The second curve shows the neutrino fluxes expected from the supernova rem- 
nant IC443 if the 7-rays detected by EGRET ^\ are a \\ Q f hadronic origin . The 
neutrino flux does not reach very high energy because the EGRET detection is con- 
sistent with a low maximum energy for the accelerated cosmic rays. The detected 



7-ray fluxes should be generated in hadronic interactions if the accelerated protons 
hit a dense molecular cloud (>300 cm -3 ) or are contained in a more tenuous cloud 
by magnetic fields. 

Neutrinos in these first two examples are generated in pp interactions while the 
next three models are based on photoproduction. Curve 3) shows the expected neu- 
trino luxes if the TeV 7-ray outburst of Mrk 501 is of hadronic origin. The 
neutrino flux extends to higher energy than the detected 7-rays as the latter may 
have been absorbed in 77 collisions either on source or in propagation to us. 

Curve 4) shows the minimum and maximum fluxes expected from the core region 
of 3C273 The photon density in the core region is estimated from the total 

luminosity of the source and the proton density is estimated from the accretion rate 
that could support the source luminosity. Magnetic field (estimated by equipartition) 
is sufficient for the acceleration of protons to high energy. Since the photon density 
exceeds the proton density by at least eight orders of magnitude, photoproduction is 
the major energy loss process for the high energy protons. 

Curve 5) shows the neutrino flux predicted for the jet of 3C279^^. Neutrinos are 
boosted to higher energy by the Doppler factor of the jet, which is 10 in this example. 

There are obviously many more, and newer, models, but the selection shown 
above includes examples of all possible types of astrophysical objects that are poten- 
tial neutrino sources. There are also galactic models that rely on photoproduction 
interactions in objects like micro- quasars. 



5.2. Diffuse neutrinos 

Figure |H1 shows several different diffuse astrophysical neutrino fluxes. The shaded 
area indicates the vertical and horizontal fluxes of atmospheric neutrinos. The Wax- 
man&Bahcall limit . derived from the flux of the highest energy cosmic rays in 
optically thin astrophysical objects is indicated with W&B. 

The curve labeled 1) shows the neutrinos expected from the central Galaxy in 
the assumption that all diffuse 7-rays detected by EGRET are created by cosmic 
ray interactions with matter. This is indeed a standard assumption for photons of 
energy above 100 MeV. The only questionable feature of this observation is that the 
7-ray energy spectrum is flatter than what is expected from interactions of the locally 
observed cosmic rays. This can be assigned to the existence of unresolved sources 
with spectra much flatter than those of the local cosmic rays. 

Curves 2) come from a cosmological integration of the models of Ref. The 
cosmological evolution of AGN is assumed to be close to (1 + z) A (although the form 
used is much more sophisticated). Since the AGN cores are optically thick these 
fluxes do not have to obey the W&B limit. 

Flux 3) is the isotropic AGN neutrino flux from Ref.^^, where pp interactions are 
added to the high energy photoproduction interactions. For about two decades(10 165 




- 10 18,5 eV) the flux exceeds the W&B limit by a small amount. The neutrino flux 
from pp interactions does not have to obey the limit. 

Flux 4) is the prediction of diffuse neutrinos from gamma ray bursts in the 
assumption that GRBs are sources of the ultrahigh energy cosmic rays. 

Flux 5) is a nominal cosmogenic neutrino flux as calculated in Ref. using the 
luminosity and cosmological evolution model from the W&B limit. 

Finally, for comparisons with diffuse astrophysical neutrinos we show the neu- 
trino flux (label 6)) that is needed by the Z-burst model to become the production 
mechanism for UHECR^l. 



6. Cosmogenic neutrinos and cosmological evolution of the astrophysical 
sources 

Cosmogenic neutrinos are generated by photoproductio n inter actions of high en- 
ergy cosmic rays, mostly with the microwave background To illustrate the 
strong cosmological link we shall build a model of the production of cosmogenic neu- 
trinos, that includes the following simplifying assumptions following Ref. **B\ 

• The microwave background is the only target for cosmic ray interactions. 

• The Universe is matter dominated (Qm = 1). 

• Neutrino production is very fast on cosmological scale. 

• Cosmic ray sources evolve as (1 + z) m forever. 



The simplicity of the model changes the result of the calculation by less than a factor 
of two. 

The neutrino yield Y from protons of energy E p then scales with the redshift as 
E U 4£-(E V , E p , z) = E v ^-(q 2 E u ,qE p ) , where Y is the yield in the contemporary 
universe and q = 1 + z. To obtain the flux of atmospheric neutrinos one can integrate 
over E p and lng and obtains 

where the yield Y 0l depends on the integral slope of the cosmic ray injection spectrum 
7- 

The interesting result here is the g( m+7 ~~2) dependence of the flux, when the result 
is expressed in an integral over lng. For (771 + 7) l ess than 1.5 the contribution to the 
cosmogenic neutrino flux decreases with redshift. If (771 + 7) = 1.5 the contribution 
of all cosmological epochs is exactly the same. We assumed Qm — to show this flat 
contribution. In a A dominated Universe a slight curvature would appear. For values 
higher than 1.5 earlier cosmological epochs dominate the flux as shown in Fig. El 




Figure 9: Cosmogenic neutrino fluxes from the simplified model described above. The thin lines 
show the contributions to the total flux from redshifts z = 0, 1, 2, ... 10, from right to left. The 
points show the total flux. 

For higher values of [m + 7) the contribution of earlier cosmological epochs would 
be so strong that they would dominate even when a more reasonable source evolution 
model (with a cut-off) is used. The total flux of cosmogenic neutrinos also significantly 
increases. 

This is directly related to the fits of the extragalactic cosmic ray spectrum. There 
are generally two types of fits: 

Scenario 1: A steep cosmic ray injection spectrum (7 > 1.5) after propagation fits 
well the measured spectrum down to 10 18 eV. The second knee*^ of the spectrum, 
a dip in the spectral shape when multiplied by E 3 , is due to the BH pair production 



by protons as predicted in^. The shape of the propagated spectrum is such that no 
cosmological evolution of the cosmic ray sources is possible. 

Scenario 2: The injection spectrum is flat (7 = 1) as expected from shock acceleration 
models. Since a flat injection spectrum can not fit the measured cosmic ray spectrum, 
a cosmological evolution of the cosmic ray sources of order (1 + z) 3 is needed. Even 
then the galactic cosmic ray spectrum has to extend up to 10 19 eV. The second knee 
is formed at the intersection of the galactic and extragalactic cosmic ray spectra. 

Both scenarios have been discussed quite actively in the literature during the re- 
cent couple of years. Each one has its supporters, although for now this is only a 
matter of preference. Scenario 1 has been presented in 1 1 1 ' and discussed in other 
recent publications. Scenario 2 was developed in ^3 and supported in and else- 
where. One difference between these two scenarios is the flux of cosmogenic neutrinos 
that they generate in the assumption of uniform and heterogeneous distribution of 
cosmic ray sources. In the case of Scenario 1 the sum 7 + m is 1.7, while in Scenario 2 
the sum is 4. Scenario 2 will thus generate much larger flux of cosmogenic neutrinos 
than Scenario 1. 
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Figure 10: Cosmogenic neutrino fluxes generated by the two scenarios for the detected flux of 
UHECR. The 7, m values are indicated in the graph. 

Figure ITU1 shows the cosmogenic neutrino spectra from these two scenarios for the 
cosmic ray emissivity of 4.5 xlO 45 ergs/Mpc 3 /yr that was derived by Waxman^. 
Cosmological model with Q M = 0.3 and = 0.7 was used for this calculation. The 
actual difference between the models could be a bit smaller since Scenario 1 requires 
higher emissivity. 

There are also other difference between the two scenarios. The main one is that 
the transition between galactic and extragalactic cosmic rays is earlier (at smaller 



energy) by one order of magnitude in Scenario 1. Since we presume that galactic 
cosmic rays at that energy are all heavy nuclei (iron) they produce different energy 
dependence of the cosmic ray chemical composition. The Auger Southern Observatory 
will hopefully measure the composition with smaller ambiguity than exists today, it 
will help to solve the problem. So would the detection of the cosmogenic neutrino 
fluxes by the current planned and constructed neutrino detectors. ANITA, IceCube, 
Mediterranean km 3 , RICE) may confirm the existence of GZK neutrinos. A next 
generation of experiments (EUSO, OWL, SalSA, X-RICE) is being planned which 
would provide sufficient statistics (10-100 GZK events per yr) to complement and 
expand the AUGER observations. Successful completion of one such experiment 
would be an important step toward understanding the sources of the highest energy 
particles in the Universe. 
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